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BARroM CLOUD EVOLUTION AND STRIATION FORMATION IN THE 
MAGNETOSPHERIC RELEASE ON SEPTEMBER 21, 1971 

By David Adamson and Clifford L, Fricke 
Langley Research Center 

SUMMARY 

The joint NASA-Max Planck Institute Barium Ion Cloud (BIC) Experiment on 
September 21, 1971, involved the release of 1.7 kg of neutral barium at an altitude of 
31 500 km at a latitude of 6.93° N, and a longitude of 74.40° W. A theoretical model 
describing the barium neutral cloud expansion and the ion cloud formation is developed. 
The mechanism of formation of the striational features observed in the release is also 
discussed. Two candidate instabilities, which may contribute to striation formation, are 
examined. The drift instability stemming from the outwardly directed drag force exerted 
on the ions by the butstrearaing neutrals is rejected on the grounds that the ion density is 
too low during the collision-dominated phase of the cloud expansion to support this kind of 
instability. The joint action of Rayleigh-Taylor and flute instabilities plausibly accounts 
for the observed striational structure. This same mechanism may well be operative at 
times of sudden injection of plasma into the inner magnetosphere during geomagnetic 
storms and may thus contribute to the formation of field-alined inhomogeneities which 
serve as whistler ducts. 


INTRODUCTION 

In the presence of a magnetic field plasmas tend to form field-alined inhomogeneities 
or striations. This tendency manifests itself in a variety of ways throughout nature. Per- 
haps the most spectacular example is provided by the striated or rayed aurora. Another 
instance is the existence in Earth’s magnetosphere of field-alined tubes of enhanced elec- 
tron density which serve as guides to the whistler mode of wave propagation. On a larger 
scale, visual evidence is noted in the filamentary structures of the solar corona and in the 
prominences which occur within it (fig. 1). On an even grander scale, the filamentary 
structures in gaseous nebulae may plausibly be argued to be manifestations of the same 
phenomena. The fact that in all of these phenomena the striations are field-alined does 
not imply a common underlying mechanism of formation. It is simply a consequence of 
the fact that in tenuous plasmas the conductivity along the magnetic field lines is substan- 
ially greater than the conductivity transverse thereto, and the same electrodynamic forces 
are therefore operative at all points along the same magnetic field line. 
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Figure 1.- Filamentary structure of solar corona. 

Photograph of solar eclipse on March 7, 1970, 
taken and processed by Leonard M, Weinstein 
of NASA LRC and Sheldon M. Smith of NASA 
ARC. 

In view of the widespread occurrence of filamentary structures in naturally occur- 
ring plasmas, it is hardly surprising that striations are invariably observed in barium 
ion cloud releases in the ionosphere. In certain of the releases at high latitude (ref, 1) 
there are strong indications that some of the striational features result from the barium 
ions acting as tracers in delineating magnetic field inhomogeneities preexisting the 
release. In other cases the regularity of the striational features coupled with their pro- 
pensity to form in the proximity of the trailing edge of the cloud leads to the speculation 
that their formation stems from some inherent instability triggered by interaction with 
the ambient medium. A number of theoretical investigations into the nature of the insta- 
bilities contributing to striation formation in barium clouds released at ionospheric lati- 
tudes have recently been made (refs. 2 to 6). The theories advanced for the explanation 
of these instabilities in ion clouds in the ionosphere all emphasize the importance of the 
role played by collisions between the drifting ions and the neutral component of the ambient 
medium. 

Two experiments involving barium releases deep within the magnetosphere have been 
made. In the first experiment (ref. 7) about 100 grams of barium were released at an 
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altitude of about 11.5 Earth radii (about 75 000 km) from the European Space Research 
Organization (ESRO) Highly Eccentric Orbiting Satellite (HEOS I) on March 18, 1969. One 
or two faint striational features were observed. The appearance of these striations in the 
almost perfect vacuum existing at such a long distance from Earth was attributed to flute 
instabilities operative during the initial phase of cloud expansion (ref. 8). The second 
experiment made on September 21, 1971, on which this paper is based, involved the release 
of 1.7 kg of barium at an altitude of 31 500 km (ref. 9). Considerable striational structure, 
consisting of about a dozen prominent striational features, evolved within the first 2 min- 
utes and remained visible for about 30 minutes. Since the purpose of this paper is to inves- 
tigate the nature of the instabilities which contribute to striation formation, consideration 
is limited to the early phase of cloud expansion. 

In the section entitled ’’Observational Data" a brief description is given of the obser- 
vational data bearing on striation formation and structure. In the section entitled "Cloud 
Model" a theoretical model, which purports to describe the early time evolution of the 
barium neutral cloud and the formation of the ion cloud, is developed. Subsequent sections 
are devoted to a discussion of drift instability and Rayleigh-Taylor and flute instabilities 
and the role they may play in striation formation. 

SYMBOLS 

a classical radius of electron 

e 

b unit vector in direction of magnetic field 

B magnetic-field vector 

strength of magnetic field outside ion cloud 
c velocity of light 

C constant appearing in expression for rate of ionization of neutral cloud 

d offset distance of line of sight from cloud center 

e electron charge 

E electric field vector 

f oscillator strength 
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F flux of radiant energy 

d^r 

g acceleration of cloud boundary, — S- 

dt^ 

G(v,t) function defined in equation (14) 

Iq differential flux of radiant energy (flux per unit frequency range) 

j electric current vector 

k wave number of boundary perturbation; Boltzmann constant 

k’ absorption coefficient per vmit volume 

K constant 

I length of cloud profile arc 

L length of ion cloud measured along cloud axis 

m electron mass 

n number density of neutral atoms 

number density of ions 

N(t) total number of neutral barium atoms as a function of t 

^N(t) number of neutral atoms having radial velocities in range v to v + Av at 

time t 

Nj(t) total number of barium ions as a function of t 

p plasma pressure 

q mode of fluting of boundary perturbation 

r radial distance from neutral cloud center; distance from ion cloud axis 
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equatorial radius of unperturbed ion cloud boundary 
R radius of curvature of magnetic field line 

R radius of curvature vector 

s distance measured along boundary in direction of magnetic field 

t time elapsed from instant of cloud release 

u velocity component of neutral barium atoms in direction of solar irradiation 

v radial velocity of neutral atoms 


Vj^ mean radial velocity of neutral atoms 

v,p radial velocity spread of neutral atoms 

Vth average thermal velocity of ions 

X axis pointing toward source of irradiation 

z distance along cloud axis 

a exponential growth of instabilities (absence of thermal quenching) 



6t) amplitude of boundary perturbation 
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K 


X 

= cy 
% 


V 


constant appearing in exponential damping factor associated with thermal 
damping of instability 

wavelength of spectral line 

circumferential wavelength of bovmdary perturbation of ion cloud 

magnetic permeability 
frequency of spectral line 





P 
a 

T 

^d 

Subscripts: 

o values at unperturbed boundary of ion cloud 

e values at instant of pressure equalization 

max maximum 

In the stability analysis, tilde over a quantity denotes its unperturbed value. 


plasma density 
variable of integration 
variable of integration 

time for ion to move between crests and troughs of fluted ion cloud boundary 
azimuthal angle 

brightness parameter ^particles per cm^ column along line of sight] 


6 



OBSERVATIONAL DATA 


On September 21, 1971, 1.7 kg of barium were released at an altitude of about 
31 500 km, a latitude of 6.93‘^ N. and a longitude of 74.40° W. The barium was produced 
by a thermite-type chemical reaction and vented into space as a neutral cloud. Under 
solar irradiation the neutral barium was almost completely ionized within the first min- 
ute. The ions, as they formed, spiraled upwards and downwards along Earth's magnetic 
field lines, the ion cloud as a whole elongating therefore along the magnetic field direction. 
Figure 2(a) is a photograph of the release taken 2 minutes 27 seconds after release with a 
Baker-Nunn camera located at Mount Hopkins, Arizona. The spherical neutral cloud con- 
sists of strontium, which was present in the payload as an impurity (the barium having 
been completely ionized at this epoch). Note that the neutral cloud exhibits the shell-like 
structure invariably accompanying the explosive expansion of neutral clouds in the hard 
vacuum of space existing at altitudes in excess of 500 km (ref. 10). The bar through the 
center of the neutral cloud comprises the main core of the barium ion cloud. The fact 
that at early times the ion cloud remains centered on the neutral cloud, which moves with 
the velocity of the rocket, results from the buildup of polarization charges within the ion 
cloud interior (ref. 11). Another feature in this photograph worthy of note is the existence 
of striational features on both sides of the main core. The fuzzy appearance of the stri- 
ations at these early times is attributed in part from limitations of resolution of the cam- 
era system used and in part from obscuration of the thin filaments by the neutral cloud in 
which they are engulfed. The photograph in figure 2(b) was taken 19 minutes 8 seconds 
after release with a 70-mm f/1 image intensified optical system located at Cerro Morado, 
Chile. At this time considerable striational structure is discernible. The regularity of 
the spacing of those striations contiguous to the main core precludes their being interpreted 
as delineations of naturally occurring inhomogeneities in the magnetic field preexisting 

I 

the release. 

In figure 3 the angular separations from the main core of all identifiable striational 
features are presented as functions of time for the period 16 minutes following release. 

This figure incorporates data obtained by using an image intensified optical system located 
at White Sands, New Mexico. For each photographic frame the coordinates of one point on 
each of the discernible striational features were read on a Telereader and automatically 
recorded on punched cards. The angular distances between the main core and the stri- 
ations were computed by using two known stars as a reference and plotted as in figure 3 
with a computer plotting system. The main core lies along the abscissa axis and at any 
epoch a section perpendicular to this axis defines the apparent angular displacement in 
degrees of the striational features relative to the main core, as viewed from White Sands, 
New Mexico. Although no reliable data exist on striation separation during the first 
4 minutes as a result of obscuration of the striational features by the overlying neutral 
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(a) Photograph taken at Mount Hopkins, Arizona (2 min 27 sec after release). 
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(b) Photograph taken at Cerro Morado, Chile (19 min 8 sec after release). 
Figure 2.- Photographs of barium cloud and component striations. 
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cloud (thou^ a 200 A interference filter was used, the closeness of the spectral lines of 
ionized barium, 4554 A, and neutral strontium, 4607 A, precluded their separation), the 
tendency of the striations to converge on the release point can be noted, thereby indicating 
that the striations did form, through some intrinsic instability, during the transient phase 
of cloud expansion. This conclusion is further reinforced by data obtained with a highly 
sensitive image-intensified television camera system (ref. 12), 

CLOUD MODEL 

The purpose of the present section is to derive expressions for the density distribu- 
tions of particles comprising both the neutral and the ionized barium clouds as functions 
of position and time. Such data are needed as input to the subsequent stability 
considerations. 

The following assumptions underlying the present analysis are listed together with 
some pertinent comments pertaining thereto: 

(1) The barium atoms comprising the neutral cloud are released simultaneously and 
move radially outward without collisions. From a prior experiment (ref. 13) the mean 
velocity of radial expansion v^ is known to be approximately equal to 1200 m/sec and the 
thermal velocity spread v,p about the mean is approximately equal to 150 m/sec. (As a 
result of the mean velocity being so much larger than the velocity spread, the neutral cloud 
exhibits the shell-like structure so evident in fig, 2(a).) By virtue of the almost explosive 
expansion of the barium neutral cloud into the hard vacuum of space, the collisions cease 
to play a role after the first second or two and can therefore be ignored. 

(2) The ionization process is dominantly controlled by absorption at a single wave- 
length. The ionization process is in fact a complex one involving a number of transitional 
steps (ref. 14), As in all such chain reactions the slowest one controls the rate of ioniza- 
tion. Therefore, for the purposes of the present calculation, ionization may be regarded 
as a single-step process resulting from absorption at the wavelength of the pacing 
transition. 

(3) The spectral line associated with the pacing transition is rectangular in form and 
has a frequency width y in hertz. 

(4) The ions, once formed, spiral along the magnetic field lines with a velocity com- 
ponent along the field line equal to that of the parent neutral atom. During ionization an 
excess of energy is imparted to the electrons, and this is subsequently shared with the 
ions; however, as a consequence of the large mass disparity, the effect on the ions is 
small. In the barium release on September 21, 1971, the rate of cloud elongation built up 
slowly to an unexpectedly high value over a period of several minutes (ref. 9). However, 
since the transient phase of cloud expansion, which is solely of interest in the present 
analysis, lasted less than a minute, this slow energization of the ions need not be considered. 
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(5) The ion gyroradius is small when compared with the cloud diameter. If the rad- 
ially outward velocity of the ions as they are formed is assumed to be the same as that of 
the parent atom, that is, approximately 1 km/sec, the ion gyroradius corresponding to a 
magnetic field strength of 150 gammas existing at the altitude of release is about 60 km. 

This is much larger than the observed cloud radius of about 8 km at later times, and on 
these grounds the validity of this assumption may, therefore, be questioned. However, 
bear in mind that the effect of the electrons, which will remain tied tightly to the magnetic 
field lines, will be to inhibit the orbital motions of the ions, thus materially reducing the 
effective ion gyroradius. This is borne out by the numerical studies of initial cloud expan- 
sion described in reference 15. 

(6) The magnetic field lines remain undistorted during ion cloud formation. At very 
early times, prior to striation formation, the magnetic field lines are probably excluded 
from the cloud interior. This is borne out by the numerical studies reported in refer- 
ence 15 in which the condition of radial symmetry is enforced. During the striation growth 
and detachment the field lines would be expected to infiltrate the cloud interior. So little is 
known about this infiltration process that it cannot be adequately incorporated into a cloud 
model. Thus in formulating a cloud model, the following two alternatives exist: (a) To 
assume complete exclusion of magnetic field lines, in which case the model would be more 
realistic at very early times; or (b) to assume the field lines to permeate the cloud inte- 
rior, in which case the model would be more realistic at later times. In the present formu- 
lation the latter course has been chosen. 

Evolution of Neutral Cloud 

In figure 4 the spherical neutral cloud is depicted as being irradiated with sunlight 
from the right. Since the atoms are assumed to be released instantaneously and move 
without collisions, those occupying the disk AA will all have velocity components in the 
x-direction in the range u to u + Au. The disk thus advances to the right with veloc- 
ity u, distending in the radial direction as it does so. The thickness of the disk at any 
instant is given by the expression Ax = t Au where t is the time elapsed from the 
instant of release. Clearly, the velocity gradient in the x-direction is uniform across the 
entire cloud and given by 

^=1 ( 1 ) 

Ax t 

Since the Doppler shift of the absorption band, and hence the ionization, is wholly 
independent of the velocity component perpendicular to x (x points toward the irradiating 
source), the problem can be conceptually simplified by treating it as a two-dimensional slab 
dilating in the x-direction, as represented in figure 5. Without loss of generality, the arbi- 
trary section XX is assumed to be at rest. The crosshatched swath superimposed in 
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Figure 4.- Expanding neutral cloud. 

figure 5 serves to define, on the frequency scale on the left, the absorptive frequency range 
as a function of cross-section location. Thus, at section XX the range of absorbed fre- 
quencies extends from to A section to the right of XX will be moving to the right 
and its absorptive waveband will therefore be Doppler shifted to lower frequencies and 
vice versa; 



Solar 

radiation 


Figure 5.- Dilating slab. 
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The Doppler shift as a function of velocity is given by 


At^ _ 1 
Au ^ 


(2) 


From equations (1) and (2), where X is the wavelength, the expression for the slope 
of the absorptive waveband in figure 5 is 


^ = 1 (3) 

Ax Xt 

The radiation, in the frequency range dt» centered on u^, incident on the plane XX 
is subjected to prior absorption only within the slab of thickness d (fig. 5). By setting 

o = Py ° and using equation (3) 

d = + CT^yXt (4) 

Denoting the unattenuated differential flux (flux per unit frequency range) of radiation, 
assumed constant over the width of the spectral line, by 1^, the radiation incident on the 
plane XX in the frequency range dv is given by 


dF = dv 


(5) 


where k’ is the absorption coefficient per unit volume. 

By assumption (3), k' is constant over the entire bandwidth y and equal to (see 
ref. 16) 


where 

n 


k’ = 


2nfe2 

mcy 


number density 


( 6 ) 


f oscillator strength 

e electron charge 


m mass of the electron 


c velocity of light 

y spectral line width 
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By substituting equations (4) and (6) in equation (5) and integrating over the entire 
line width, the expression for the integrated flux at section XX is 


p+ 1/ 2 

Integrated flux = L \ exp - 

° J-1/2 


2nfe Xtfl . J 


2nfe2xt \ me . 


A comment is called for at this point regarding the derivation of equation (7). Note 
that implicit in equation (5) is the assumption that the density is uniform throughout the 
absorptive layer, whereas in the actual case of the barium cloud density gradients do exist. 
Furthermore, the outermost layers of the expanding spherical cloud are not screened by 
overlying layers as is assumed in the foregoing analysis. For equation (7) to be applica- 
ble to the neutral barium cloud, the thickness of the absorbing layer must be small com- 
pared with the thickness v,j,t of the spherical shell of the cloud within which the bulk of 
the barium resides. 

The total thickness of the absorbing layer contributing to flux attenuation at a given 
section, as obtained from equation (4) by setting a = I, is equal to yXt. Thus the ratio 
of thickness of absorptive layer to shell thickness is given by yX/vrp. 

^ 2p2 

By substituting the expression — — %- (see ref. 16) for spectral line width y and 

3mc3 2 

making use of the fact that the classical radius of the electron a is equal to 

® mc^ 

Thickness of absorptive layer _ /o\ 

Shell thickness SXvijr 

By assuming a Vrp of 150 m/ sec and a X of 5000 A, this ratio is found to be less 
than 0.05, which is small enough to justify using equation (7) throughout the following 
analysis. 

The rate at which neutral atoms are being ionized, and hence lost to the neutral popu- 
lation, is proportional to the radiation flux and the neutral particle density n. Thus 

= -CFn foj 


By substituting equation (7) 


2fe^Xt 


1 - expl 


2nfe^Xt 
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where 


fi = Cl^y 


Denote by AN the number of neutral atoms having radial velocities in the range v 
to V + Av. At the instant t these atoms occupy a spherical shell having an inner radius 
tv and a thickness t Av. Hence from equation (10) the instantaneous rate of ionization 
per unit volume is given at all points within the spherical shell by the expression 


■func 


2fe2xt 


1 - exp 


' o ^ 

-fe^X AN 

,27rmct^^Avy 


from whence it follows that by multiplying by the volume of the shell 


d(AN) _ -nmc 
dt 2fe2\t 


1 - exp 


-fe'^^X AN 




27rmct2v2Av^ 


4jrt^v^Av 


= -Mt 


2 2irmcv Av 


fe^x 


1 - expl 


-fe^X AN 
27rmct2v2Av> 


( 11 ) 


From an analysis of neutral cloud expansion in other high -altitude neutral barium 
cloud releases (ref, 13), the initial velocity distribution has been found to be of the form 



AN = K jexp 

Thus, the initial total number of neutral atoms N(0) is given by 

.21 


( 12 ) 


N(0) = K r 

Jo 


exp 


V - V, 


v^dv 




where 



and 


P(^) = 


— 

,1/2 Jo 


e-x dx 
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Since ^ is approximately equal to 8, then to a sufficient accuracy 
N(0) 


K = 




By introducing the constant into equation (12), it follows that at t = 0, the total num- 
ber of neutral atoms in the range v to v + Av is given by 


AN(0) = 


N(0) 


1/2 /.2 ^ l\ 3 

” ^ jjvx 


exp 


V-V^x3 


V^Av 


The function 


G(v,t) = 

2rmcv^Av 


(13) 


(14) 


is introduced and substituted into equations (11) and (13). Thus 


.^2 

dt ^ 


l-expl-^ 


'J 


and 


G(v,o) = 


2w^/^mc 


exp 


V - V, 


J 


(15) 


(16) 


Equations (15) and (16) determine G(v,t). The total number of neutral atoms is obtained 
by summing over all subsets 

N(t) = y AN(t) = y 27?mcv^Av 
^ ^ fe2x 


In the integral limit 


N(t) = ^ r 


fe^X 

The total number of barium ions 


(17) 


N.(t) = N(0) - C G(v,t)v^dv 

fe^x *^0 


(18) 
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As mentioned earlier from an analysis of the data obtained on a previous high- 
altitude release (ref. 13), the mean radial velocity of the neutral barium atoms v^ and 
their velocity spread Vrj. have been found to be approximately equal to 1200 m/ sec and 
150 m/sec, respectively. In addition, the amount of barium released on September 21, 
1971, has been estimated at 1.7 kg (n( 0) = 7.45 x 1024). This estimate is based on an 
analysis of photometric data (ref. 17). The parameter p. appearing in equation (15) has 
been determined from earlier investigations (ref. 13) to be equal to 0.05. The remaining 
parameter fX depends on the nature of the pacing transition involved in the ionization 
process, which remains imidentified, and must be regarded as an unknown. By varying 
this parameter and matching the predicted time dependence of the total number of ions as 
calculated from equation (18) against the observed time dependence as deduced in refer- 
ence 17 from the photometric variation of ionized cloud brightness at 4554 A, it was found 
that the best fit over the initial 50 seconds was obtained by assigning fX = 1500 A. The 
excellent agreement obtained is clearly exhibited in figure 6. 



Figure 6.- Fraction of ions formed as function of time. 

At epoch t the neutral atoms populating the spherical shell of inner radius r and 

A Y* 

thickness Ar have velocities in the range v = ^ and Av = -j-. 
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Hence, from equation (14) the number of neutral atoms between r and r + Ar at 
epoch t is equal to 


2ffmci 


(fff 


fe^X 


It ’ ) 


Thus the density of neutral atoms at r,t 


is given by 


The neutral density obtained from this equation is plotted against r for the four 
epochs, 10; 20, 30, and 40 seconds in figure 7. Note that the density peak propagates out- 
wards at a rate equal to the mean radial velocity of the neutral atoms of 1200 m/sec. Note 
also that the vertical scales are different at the four epochs. This difference was neces- 
sitated by the rapid falloff of neutral density resulting in part from the cloud expansion 
and in part from the rapid ionization. 

Variation of Density Within Ion Cloud 

For simplicity only the ion distribution in the median plane, drawn throu^ the cloud 
center perpendicular to the field lines, will be considered. The density of neutral atoms 
at instant r in the proximity of the point P a distance r from the origin (fig. 8) is 
given by equation (19). By noting that the rate of formation of ions is equal to the rate of 
decrease of the neutral population, the following expression for the number of ions depos- 
ited within unit volume centered on P in the interval dr is obtained from equations (19) 
and (10) 



At the instant of deposition the velocity gradient across the neutral cloud is equal to 
1 /t (see eq. (1)) and by virtue of assumption (4), that the velocity component of the ions 
in the direction of the magnetic field lines is the same as that of the parent neutral atom, 
it follows that the volume element in which the ions are deposited at instant t will there- 
after continue to elongate at a rate of 1 /t. Thus at epoch t the length of the unit volume 
(unit cross section and xmit length) will have become 

1 + i (t - t) = 1 

T ' T 


18 



Radius, r, km 


Radius, r, km 


Figure 7.- Radial density distribution of neutral barium atoms as function of time. 









Figure 8. - Elongation of volume element with time. 


and thus the ion density contributed at instant t by those ions formed during the inter- 
val T to T + At will be equal to 


dn. = 

^ 2fe2xri 


1 - exp 



( 21 ) 


By integration the following expression for the density of ions in the median section 
at radius r at instant t is obtained 


Curves showing the variation of ion density with time at radii of 5, 10, and 15 km in 
the median section are plotted in figure 9. Note that the density at a given point remains 
small until the high -density portion of the neutral gas shell encroaches upon it. This is 
followed by a period of rapid buildup of ions as the dense shell passes over it, depositing 
ions as it does so. Once the dense portion of the neutral gas shell has passed by, there is 
a slow falloff in ion density as the ions disperse along the field lines. Figure 10 shows 
the radial ion distribution within the median section at epochs of 10, 20, and 40 seconds. 
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Density of barium ions/cm 





Figure 10. - Radial distribution of barium ions in median cross section 

at selected epochs. 


Brightness Variation Across Median Section of Ion Cloud 

The brightness of an optically thin cloud is proportional to the number of ions within 
a column of unit cross-sectional area. By assuming the cloud to be viewed at instant t 
along a line of sight offset by a distance d from the cloud center (fig. 11), the brightness 
is given by 


K 



(23) 


where n.(r,t) is given by equation (22). Note that the integral in equation (23) defines the 
number of ions along a line perpendicular to the cloud axis. In the event the cloud is 
viewed obliquely, then the number of ions along the line of sight is equal to the integral 
(eq. (23)) multiplied by where 0 is the angle between the line of sight and the 

cloud centerline. This factor is assumed to have been absorbed in the constant K of 
equation (23). 

The expression 
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Figure 11. - Median cross section of ion cloud. 

has been calculated as a function of p for the epochs of 10, 20, and 40 seconds following 
release. The results are presented in figure 12. Note the indentation of the brightness 
profile at early times in the proximity of the cloud axis. This could not be checked obser- 
vationally because of obscuration of the ion cloud at early times by the superimposed 
neutral cloud. 


DRIFT INSTABILITY 

In the drift instability which leads to the formation of striations in ionospheric 
releases, an important role is played by the collisions taking place between drifting ions 
and background neutral particles, hence the name drift instability. There has been spec- 
ulation that the same kind of instability may be operative at much higher altitudes, where 
the background density is virtually zero, by invoking collisions between outstreaming 
neutral barium atoms and the barium ions. Therefore a brief description of the insta- 
bility and the reasons for rejecting it as a contributory cause of striation development in 
magnetospheric releases is given. 

The underlying mechanism is illustrated in figure 13 depicting a cross section of the 
field -alined ion cloud relative to which the background neutrals are streaming from left to 
right (to an observer at rest relative to the neutral population, the cloud is drifting to the 
left and hence the forward and rearward faces are as indicated). If the magnetic field is 
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Figure 12.- Profile of brightness parameter $ across median section of ion cloud. 

assumed to be out of the plane of the paper, the drag force exerted by the neutrals on the 
charged particles will result in a drifting of positive ions and electrons in the downward 
and upward directions, respectively. If it is further supposed that the forward and rear- 
ward faces are rippled, as shown in figure 13, the accumulation of charge and the associ- 
ated electric field E within the plasma will be in the direction depicted. Clearly, the 
polarization electric field causes the plasma to drift with velocity E x b/b 2, which is in 
a direction to suppress the undulation on the forward face and to amplify it on the rear- 
ward face. 

Turning now to a consideration of the barium releases in deep space, the ions as 
they form are constrained by the magnetic field and at early times there is a radial out- 
streaming of neutrals through them, as depicted in figure 14. Clearly, the conditions 
existing at the surface of the ion cloud correspond to those prevailing at the rearward 
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Figure 13.- Cross section of ion cloud released in ionosphere. 



Figure 14.- Cross section of ion cloud released in deep space. 
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face of the ionospheric cloud in the previous example (fig. 13). Conditions are, therefore, 
conducive to instability. 

By using the results of the analysis of the preceding section, a plot has been made 
of peak neutral density as a function of time (fig. 15). By adopting a collisional cross- 

a 

sectional radius of 6 A for encounters between neutral barium atoms and barium ions, the 
interval between successive encounters of an ion with neutrals is given as a function of 
time in figure 16. Note that at 1.5 seconds after release the interval between collisions 
is 2.0 seconds, and after 5 seconds the collisional interval has increased to about 75 sec- 
onds. Obviously collisions between ions and neutrals can be ignored after several sec- 
onds. However, referring to figure 6, note that too few ions are formed within the collision- 
dominated phase of the cloud expansion to support the type of instability under discussion. 



Time, sec 


Figure 15. - Maximum density of neutral cloud as function of time. 
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Figure 16.- Interval between collisions of an ion with neutrals as function of time. 

RAYLEIGH-TAYLOR AND FLUTE INSTABILITIES 

In the present section the assumption is made that the rapid formation of the ion 
cloud in the weak magnetic field existing at 5 Earth radii results in a thrusting apart of 
the magnetic field lines. Thus, whereas the cloud elongation in the direction of the mag- 
netic field proceeds without impediment, its transverse expansion is opposed by an inwardly 
directed magnetic pressure. During the resulting radial deceleration the inertial forces on 
the plasma are outwards. The situation is clearly analogous to one in which a heavy fluid 
is supported by one of zero density in a gravitational field, and the cloud boundary is there- 
fore prone to the so-called Rayleigh -Taylor instability. In addition, the configuration is 
subject to flute instabilities of the kind discussed in reference 18. Thus, suppose that the 
cloud boundary becomes fluted, as shown in figure 17, and that the magnetic field is wholly 
excluded from the cloud interior, then the integral of the magnetic field vector taken around 
a circuit extending from a to b along a crest (solid line in fig. 17) and from b to a 
along an adjacent trough (dashed line in fig. 17) is zero since no current is threading the 
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Figure 17.- Ion cloud with fluted boundary. 



circuit. It follows therefore that since the arc length along a crest exceeds that along a 
trough, the magnetic field pressure at a crest is less than that at a trough, a condition con- 
ducive to further growth of the flutes. Both of these instabilities, Rayleigh-Taylor and 
flute, are quantitatively analyzed in this section. The general approach used has similari- 
ties with that adopted in reference 8, however the manner of derivation as herein presented 
differs from that of the aforementioned reference as do the end results. 


Analysis 

By introducing cylindrical coordinates, the axisymmetric unperturbed boundary is 
given by the expression r = r(z,t). In the stability analysis consideration is limited to 
motion in the equatorial (median) plane. Denote the unperturbed equatorial radius of the 
ion cloud boundary by r^^, which is a function of time t, and denote the radial distance 
from the unperturbed boundary by tj measured positive outwards. The equation of motion 
in the radial direction is given by 


2 


dt2 
d^Tj 


(Txb) 


V 


where 


d^r 


g = -■ 


m 


dt2 


In the magnetohydrodynamic approximation, which is used throughout 


(24) 


1 


— V X B 


Thus 

f XB = -Abx (vxb) = +X(b . v)b (25) 

If b denotes the unit vector in the direction of the magnetic field, then the operator 
b • denotes differentiation along the magnetic field lines. Bearing in mind that in 

the proximity of the median plane the magnitude of the magnetic field is constant to a first 
order, it follows from figure 18 that 



i = Bf =-B 


2R 

R 


where R is the radius of curvature vector of the magnetic field lines. 


(26) 
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Figure 18.- Cross section through ion cloud axis. 

As stated previously, at very early times prior to striation formation exclusion of 
magnetic field lines from the ion cloud interior would be expected. If this is assumed to 
be true, then equation (24) becomes 

(271 


External to the barium cloud the ambient plasma pressure and density are to all 
intents and purposes zero. Bearing this in mind and substituting equations (25) and (26) 
into equation (24) 


0 = - 4 :^ 




where denotes the magnetic field exterior to the cloud. 

Equation (28) expresses the balance of forces acting on the magnetic tube seg- 

ment aa^bb’ in figure 18. The inward force ^ results from the tensions acting at 

u R 

° /b 2 \ 

II . • 1 -1 T • -1 1 1 ^ X 1 ^ f e 1 G G i 


the ends aa^ and bb*, and the outward force results from the compressive 

pressure exerted by contiguous magnetic tubes across the lateral faces ab and a'b . 
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When the values of the unperturbed parameters are denoted by the addition of a tilde 
over the appropriate symbol, the unperturbed pressure is denoted p(? 7 ,t), and so forth. 
With this notation equations (27) and (28) become 


P = Pq + pgr? 


(29) 


where the subscript o corresponds to the location of the unperturbed boundary, that is, 
7] = o, and 





= 0 


The pressure balance at the unperturbed boundary further requires 


(30) 


P 


o 



(31) 


Now assume a perturbation of the boundary having the form of magnetic field-alined 
ripples as shown in figure 17. Thus, the perturbed boundary in the equatorial plane is 
given by the expression 


T] = 


(32) 


where 0 is the azimuthal angle. 


Figure 19 depicts the boundary distortion in the equatorial plane 

2j7rm 

The circumferential wavelength of the fluting is given by X = . 

^ q 

wave number is given by 


for the mode q = 8. 
The corresponding 


^c ^m 

The perturbation to a field quantity 
has the form 


(33) 

P associated with the boundary perturbation 


6(P) oc 
6(P) oc 


(outside the boundary)] 
(inside the boundary) [ 


(34) 


where the sign of the exponent defining the radial dependence of the perturbation is so 
chosen to ensure an exponential falloff in the perturbation away from the boundary in 
either direction. Specifically 
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Figure 19.- Equatorial boundary perturbation (q = 8). 


P = Pq + Pg^ + 6pe^^e^^^ 'I 



(35) 


The radial motion of a plasma element contiguous to, and inside the boundary, will 
necessarily be the same as that of the boundary* By substituting equations (32) and (35) 
into equation (27) 

+ pgTj 4- 6pe*^%^^^)+ pg 

Thus 

pdj^ = -k6p (36) 

dt^ 

By substituting equation (35) into equation (28) and subtracting equation (30) to first 

order 
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Figure 20.- Longitudinal cross section of ion cloud. 
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From Longmire (ref. 19) 

5B„ 


where 


6B, 


61 


' I 


perturbation in magnetic field at crest in equatorial plane 


amount by which length of crest profile exceeds undeformed profile 


By assuming the equation of the undeformed profile to be 






\2 

j dz 


its profile arc length measured between points a and b (fig. 20) is given by 
I 

For the crest profile 

r + dr = (1 + e)-|i^l - cos 

and its arc length between points a and b is therefore 


I + 




( 


dz^ 


V \dzJ 


:dz 


( 41 ) 


(42) 


(43) 


(44) 


(45) 


From equations (43) and (45) and by assuming the cloud profile to be reasonably slender 


dl 


L / \2 
J dz 




By substituting equation (42) 


rL 

dl ~e' “ 




dz 
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Hence 



Moreover 


Ro 



( 46 ) 


(47) 


By substituting equations (41), (46), and (47) into equation (40) and bearing in mind 


that 5ri - er 


m 


2 


P^-kPg + lH— £ 


m 


dt" 


2 IIq l 2 


6/? 


,(Vt 


By assuming an exponential growth of the instability 6rj ^ e and substituting 
expression (33) for k, the dispersion relation becomes 


a = \/qi 


-g_ + 


2 B, 


(48) 


This relation implies that the exponential growth rate increases without bound as q 
increases. This, however, is physically unrealistic. Remember that the ions can move 
to and fro within the field free cloud interior, and therefore their thermal motion will tend 
to erode the crests and fill in the hollows of the fluted boundary. The time taken by an 
average ion to move from a crest to a hollow as a result of thermal motion is clearly 


TTr, 


= 


m 




(49) 


In the absence of instabilities a subsidence of flutes resulting from thermal agitation 

-fct /Td 

as given by the expression e where k is a numerical factor might be expected. 

The net growth rate is thus given by 

e«t ^ gOitg-'^V^d 


Thus 



(50) 
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For sufficiently small wavelengths (large enough q) the exponent a becomes neg- 
ative and the flutes are eroded more rapidly than they can amplify through inherent insta- 
bility. Thus, in this case, thermal agitation of the ions quenches the instability. 


For convenience introduce 


^ = 




2 M-o pl2 


(51) 


Hence 





(52) 


The mode of deformation subject to greatest growth rate is obtained from the equation 


^ = iji. = 0 

dq 2 iq 


q = 




(53) 


and the corresponding maximum growth rate is given by 


a 


max 


7T)3 r 


m 


4kv, 


th 


(54) 


This analysis is based on the premise that there is complete exclusion of magnetic 
field from the cloud interior. This is at variance with the assumption invoked in the analy- 
sis of cloud evolution given in a previous section. The true answer lies somewhere between. 
Numerical simulations of cloud release and ionization described in reference 15 support the 
contention that there is complete exclusion from the inner core of the ion cloud; however, 
in these numerical studies the condition of radial symmetry was enforced and no allowance 
was made for three-dimensional effects. The presence of a magnetic field in the cloud 
interior would not be expected to influence markedly the contribution of the Rayleigh- 
Taylor-type of instability; however, it would increase the growth time of the flute-type 
instability. In addition, the reduced constraint on cloud expansion would reduce the effec- 
tive acceleration g and an interior field would inhibit thermal motion of the ions, thus 
modifying the damping factor. The net result of all of these competing effects is difficult 
to predict. 
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The buildup of electron charges on the faces of the surface flutes gives rise to polari- 
zation electric fields within the cloud interior. These electric fields propagate downwards 
along the magnetic field lines at Alfven velocity and are shorted out at ionospheric levels 
where the Pederson conductivity is large, foi the context of ionospheric releases these 
effects result in a partial quenching of the instabilities (ref. 4). For the barium release 
at 31 500 km there is a lapse of about 25 seconds before the polarization field manifests 
itself at ionospheric levels. Bearing in mind the rapid growth of the instabilities indicated 
by both observation and results of numerical studies given in the following subsection, the 
time of about 50 seconds before such ionospheric perturbations can rebound to the cloud is 
too long to exert any significant influence. 


Applications to Barium Ion Cloud Release at 31 500 km 

If complete exclusion of magnetic field lines from the cloud core is assumed, then 
to achieve a pressure balance at the core boundary 


where 



n„kT 

e 


(55) 


B ambient magnetic field strength (l50 gamma = 15 x 10"^ tesla 

ILC^ = 47t X lO""^ H/m 


n equilibrium ion density 


k Boltzmann constant 

T initial ion cloud temperature (4000 K) 

Thus, from equation (55) an equilibrium ion density of n^ = 1.63 x 10^^ ions/m^ is 
obtained. 

The initial rate of increase of cloud length between half -brightness points is given 
in reference 9 as 1400 m/ sec. The overall rate of elongation might be expected to be 
about twice as great. Thus, by using subscript e to denote values at the instant of pres- 
sure equalization 


L.(meters) = 2800t 
0 0 


(56) 
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The initial rate of radial expansion of the ions may plausibly be taken equal to the 
mean radial velocity of the parent neutral atoms (1200 m/sec). This radial velocity drops 
to zero as equilibrium is approached. Thus, the rate of radial expansion averaged over 
the entire deceleration period is 600 m/ sec, and 

r„ (meters) = 600t (57) 

0 C 

The total barium yield is estimated to have been 1.7 kg(N(0) = 7.45 x 10^4 barium 
atoms). If the ion cloud is assumed to have the form of a prolate spheroid and the number 
of ions to be given by the simple relation Nj(t) = N(0)(l - then at the instant of 

pressure equalization 

Ivr^^LgOg = ^1 - e x 10^^ (58) 

By substituting the value n = L63 x 10^^ ions/ m3 and equations (56) and (57) into 
equation (58), the following equation is obtained 

te^ 4 

TOTOSP - 2.16 X 10 

1 -e ® 

From this equation a deceleration time of 25 seconds is computed. Hence, the cloud 
length at midpoint of deceleration phase is L = 2800 x 12.5 = 35 000 m. 

By takin g v^ = 1200 m/sec and r^ = 600 x 25 = 15 000 m, for the effective 
cloud deceleration 

2 

V 

g = = 48 m/ sec2 

2^6 

Thus, the cloud radius at epoch t is given by 
r = v^t-lgt2 

and the cloud radius at midpoint of deceleration phase is 

r^ = 11 200 m 
m 

By introducing these numerical values into equation (51) 

/3 = 14.1 X 10"^ sec'^ 
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The mean ion velocity Vj.|^, corresponding to an initial ion cloud temperature of 4000 K, 
is equal to 900 m/sec. By using these values in equation (53), the equation obtained for 
the harmonic mode most susceptible to instability is 



(59) 


The exponential growth time associated with this most unstable mode is from equa- 
tion (54) equal to 7.3 k seconds. In figure 2(b) about 10 prominent striational features can 
be discerned. Setting q = 10 in equation (59) yields k = 0.74 and an exponential growth 
time of about 5.3 seconds. 

From equation (52) the variation in exponential growth time a with harmonic 
mode of azimuthal perturbation q is determined. The results are plotted in figure 21 
for K = 0.74. 



Mode of azimuthal perturbation, q 


Figure 21.- Exponential growth time as function of mode of azimuthal perturbation. 
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CONCLUDING REMARKS 


The regularity of the striational features observed in the barium ion cloud (BIC) 
release on September 1971 precludes their being attributed to the delineation of field- 
alined inhomogeneities preexisting the release. The observational data further indicate 
that the striations did form at early times, presumably during the transient phase of cloud 
formation. 

Throughout the initial phase of cloud evolution there is a radial outstr earning of 
neutral atoms through the ions, which have been formed and are constrained by the perme- 
ating magnetic field. The conditions prevailing at the ion cloud boundary are thus quite 
similar to those existing at the rearward face of ion clouds formed at ionospheric altitudes 
and drifting with respect to the background neutral atmosphere. For the ionospheric 
release these conditions are conducive to the onset of drift instability and are believed to 
play an important part in striation formation. The possible role played by drift instability 
in formation of striations in magnet ospheric releases has been examined. The ion density 
is too low during the collision-dominated phase of cloud expansion to support an instability 
of this kind. 

The other instability investigated in the present paper, the so-called Rayleigh- Taylor 
flute instability, is believed to explain plausibly the striational structure observed in the 
BIC experiment, both with regard to number of striations formed and their formation rates. 
Essentially the same instability had been invoked by Pilipp to explain the formation of 
some tenuous striations in an earlier magnetospheric release of March 18, 1969. Although 
there appears to be agreement regarding the basic mechanism of instability, the analytical 
treatment presented herein and the dispersion relation describing the instability growth 
differ from those of Pilipp. 

The question as to whether there is exclusion of magnetic field lines from the ion 
cloud interior during its expansion remains an open one. In the present paper two contrary 
viewpoints have been adopted; thus in the treatment of ion cloud formation and the drift 
instability it has been assumed that there is no distension of the magnetic field during the 
ion deposition, whereas in the discussion of the Rayleigh- Taylor flute instability it has 
been assumed that there is complete magnetic field exclusion from the ion cloud interior. 
The true answer probably lies somewhere in between. 

The destabilizing mechanism proposed for the striation formation in the present 
experiment would presumably be operative at times of sudden injection of plasma into the 
inner magnetosphere during geomagnetic storms. This may well contribute in part to the 
formation of those field-alined inhomogeneities in the magnetosphere which serve as 
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whistler ducts. In this connection it would be of interest to determine whether the number 
of whistler ducts tends to correlate with geomagnetic storms. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., July 1, 1974. 
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